IUHET-276 
April 1994 
Revised May 1994 



NEUTRINO FLAVOR EVOLUTION 
NEAR A SUPERNOVA'S CORE 

J. Pantaleone 

Physics Department 
Indiana University, Bloomington, IN 47405 

ABSTRACT 

In supernovae and in the early universe, neutrino flavor evolution is a 
many-body phenomena. Here the equations describing the evolution of the 
density matrices in phase space are derived. Then these equations are ap- 
plied to neutrino emission from a supernova core. The effects of a 'small' 
background neutrino density on adiabatic and nonadiabatic flavor evolution 
are calculated analytically. It is found that when flavor evolution is sizeable, 
the sensitivity to the small neutrino background is enhanced. This implies 
that r-process nucleosynthesis in supernovae may not reliably probe neutrino 
masses less than about 25 eV. 
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Forward scattering can have an enormous impact on neutrino flavor evo- 
lution. Neutrino forward scattering off of an electron background is respon- 
sible for the Mikheyev-Smirnov-Wolfenstein (MSW) effect [|TJ |2| which may 
explain the solar neutrino problem (for a review, see e.g. [0]). Neutrino for- 
ward scattering off of neutrinos @, || || has been shown to produce curious, 
coherent phenomena in the early universe @ . Here I consider flavor evolution 
when the electron and neutrino densities are comparable. This is the case 
just outside the core of a supernova where nucleosynthesis of heavy elements 
is believed to occur. Since nucleosynthesis is very sensitive to the relative 
neutrino spectra, it may be possible to probe a range of neutrino masses and 
mixings |§ which are of great interest as candidates for nonbaryonic dark 
matter. 

This paper is divided into two parts. First the general equations describ- 
ing neutrino flavor evolution in phase space are derived in the density matrix 
approximation (see also @, [10], [7[ [11], |12|). All of the conceptual steps are 
given here, but the mathematical details will be published elsewhere. These 
evolution equations have many obvious applications in the early universe and 
in supernovae. In particular, in the second half of the paper they are solved 
to find the corrections to level crossing from a 'small' neutrino background. 

1 The evolution equations. 

One starts be making the physically reasonable assumption that each neu- 
trino is described by its own wavefunction (the Hartree approximation). That 
is, each neutrino obeys a Dirac (or Majorana) equation with a potential. The 
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potential for a particular neutrino is calculated by adding up all possible two- 
body, forward scattering, weak interaction amplitudes between that neutrino 
and all the other particles present. Note that these initial two steps neglect 
Pauli blocking and so inherently assume that the neutrinos are dilute. 

Then each individual neutrino wave function is factored into a 'quickly' 
varying part (the solution to the massless wave equation) and a 'slowly' 
varying part. The 'quickly' varying part describes the propagation of a wave 
packet, and the 'slowly' varying part describes how the wave packet gradually 
separates into different wavepackets — one for each mass eigenstate. Because 
practical position 'measurements' of neutrinos are always much much coarser 
than the size of the wave packet (which is the inverse of the neutrino momen- 
tum spread), these individual equations should be spatially averaged over a 
scale intermediate between the 'quick' and 'slow' scales. This eliminates the 
interference effects from the 'quickly' varying part of the wavefunction, so 
that the 3-dimensional position and momentum vectors are then classical 
variables. There remains approximate equations for the 'slowly' varying part 
of the individual neutrino wave functions. These describe the interference 
effects resulting from the different neutrino vacuum masses. 

The potentials depend on bilinears of the fermion wavef unctions, i.e. par- 
ticle densities. With a little algebra, the equations for the individual neutrino 
wavefunctions are easily converted to equations for the bilinears 
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Pi(x,t,k) = 



V e v T 



(x,t,k) 



(1) 



Here v a is the slowly varying part of the ith neutrino's wavefunction of flavor 
a with momentum k at time t and position x (a bar over a quantity always 
denotes that it refers to an antineutrino). The individual evolution equations 

— * 

for a given x, t and k are identical. Thus I can sum the p^s and pYs over 

— * 

all neutrinos with momentum k at time t in the 'small' volume at position x 
(over which the full equations where averaged) to get the density matrices 



p(x,t,k) = ^2pi(x,t,k) , p(x,t,k) = J2pi(x,t,k) 



(2) 



The integral of p aa over all phase space is just the total number of neutrinos 
with flavor a. The coupled evolution equations for the density matrices are 



d - d ) - - 

+ k- — \ p(x, t, k) = [V(x, t, k; p, p),p(x, t, k)] 



dt 



and 



[ d - d ) _ - - -> 
M ^7 + k - ^ ( > *> k ) = t, k; p, p),p(x, t, , 



(3) 



(4) 



x dt dx 

Here the right hand sides are the commutator of the matrices. The effective 
potential matrices V and V are given in the flavor basis by 



V(x,t,k;p,p) = 
m 



(5) 



2k U 



2 



m\ 



m| 



C/ f + V2G F [N e (x, t) - N s (x, t)\ 



1 





(27T) 



and 



V(x,t,k;p,p) = 

m\ 

[^U* m\ U T + V2G F [-N e (x,t) + N s (x,t)] 

m\ 



(6) 
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The first term in the potentials describes the neutrino mixing from vacuum 
parameters, where the mf are the neutrino vacuum masses and U is the lep- 
tonic analogue of the Cabibbo-Kobayashi-Maskawa mixing matrix (see Ref. 
H and references therein for a discussion of issues relevant for 3 neutrino 
flavors). The second term in the potentials comes from neutrino forward 
scattering off of nonrelativistic electrons via the weak charged current, where 
N e (N s ) is the electron (positron) number density Jl[. Neutral current for- 
ward scattering of neutrinos off of electrons, protons and neutrons in the 
background gives a contribution proportional to the identity matrix — which 
has no effect on the flavor evolution and so is not shown. 

The final term in the V's comes from neutrinos forward scattering off 
of background neutrinos and antineutrinos via the weak neutral current |J. 
This term does contribute to the flavor evolution because the 'propagat- 
ing' and 'background' neutrinos and antineutrinos can exchange flavor. This 
term, like all tree level neutral current interactions, is invariant under under 
global flavor transformations of the form 



p' = W ] P W , p' = W T pW* 



(7) 
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where W is a unitary matrix. If p and p are in the flavor basis, and p' and p' 
are in the vacuum mass basis, then W is just U. Note that the (1 — k-q) factor 
in the neutrino-neutrino forward scattering term means that an individual 
neutrino can not forward scatter off of itself. 

These equations describe the effects on neutrino flavor evolution from 
background neutrino and electron densities. They are necessary when neu- 
trino densities are sizeable such as in supernovae and in the early universe. 
However note that nonforward neutrino scattering has been neglected in these 
evolution equations (see || [1C], |TTJ and references therein) . Thus these equa- 



tions are appropriate for supernova neutrinos outside of the neutrinosphere, 
and for cosmological neutrinos after neutrino decoupling (where additional 
CP even terms of order G 2 F are sometimes relevant, see ||). Note that these 
evolution equations do not make the two flavor approximation, and do not 
assume that phase space is isotropic. 

Under CP symmetry, p a (x,t,k) <-> p a (—x,t,—k). The wave equations 
are CP symmetric if U = U* and (N e — N g ) = 0. 

The left hand side of Eqs. (|3|-[D is just the total time derivative of the 
density matrix (since k is just the velocity of the wave packet). When all the 
quantum interference effects average out, the right hand side of these equa- 
tions vanish and so they reduce to the continuity equations in phase space, 
i.e. the Vlasov equation. The Vlasov equation generally also contains deriva- 
tives of the neutrino density with respect to momentum. Such derivatives are 
also generally present in these equations too, in principle. However changes 
in neutrino momentum come from nonforward scattering, or refraction from 
density variations, or from gravitational redshift. These effects are higher 
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order and have been consistently neglected above. 

The right hand side of Eqs. vanishes when one takes the trace of 

the matrix. The traces of p(x, t, k) and p(x, t, k) are the flavor independent 
number densities in phase space and these are conserved quantities — they 
satisfies the continuity equation in phase space. The trace of V also does 
not contribute to the right hand side of these equations. Thus in studying 
neutrino flavor evolution, I can always choose to work with a traceless p and 
p. In fact, it is easily shown that Tr[p 2 (x, t, k)] and Tr[p 2 (x, t, k)] are also 
conserved quantities. This is because dissipative effects such as nonforward 
scattering have not been included in the evolution equations, so entropy is 



conserved JL^, [TTJ] . Taking into account all of the constraints, the flavor 
content of the three-flavor density matrices depends on seven independent, 
real variables at each point in phase space. In the two-flavor approximation 
the density matrices depend on only two independent, real variables. 

It is usually true that the neutrinos are created at high densities, where 
the vacuum neutrino masses are negligible compared to the induced poten- 
tials. Then the neutrinos are approximately in flavor states. Because the neu- 
trinos are generally produced over times and positions which are distributed 
much more broadly than the scale set by the inverse of the potential, the 
flavor off diagonal terms generally average out in the summation of Eq. @. 
A density matrix diagonal in flavor space is a solution to Eqs. when 
the neutrino masses are negligible. This is typically the appropriate initial 
condition for neutrinos created at high densities. 

With the simplifying approximation that there are only two neutrino 
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flavors, the effective potentials in the flavor basis take the form: 



V(x,t,k;p,p) = 

-C20 5*26) 

d?q 



A 

4k 



/2 

+ ^-G F [N e (x,t)-Ne(x,t)} 



1 
-1 



+V2G F J - k ■ q) (-p(x,t,q)* + p(x,t,q)) 



and 



(8) 



V(x,t,k;p,p) = 

-C2Q O20 
S20 C20 

d 3 q 



A 

4k 



[2 

+ ^-G F [-N e (x,t)+N,(x,t)} 



1 
-1 



W2G F J -^(l4-g) (p(x,t,q)-p(x,t,q)*)} 



(9) 



(27T) 

here C20 = cos 29, S28 = sin 26, where 9 is the vacuum mixing angle, and 
A = ml — m\. 



2 Supernovae Nucleosynthesis. 

In a recent Letter ||, Qian et al. showed how the nucleosynthesis of heavy 
elements in type II supernovae is sensitive to neutrino masses and mixings. 
Nucleosynthesis would not occur if the v e and the or v T fluxes emitted from 
the supernova core were interchanged before reaching the nucleosynthesis 
region. Such a flavor interchange occurs when a massive neutrino propagates 
adiabatically through an MSW resonance. Qian et al. calculated that this 
occurs for neutrino mass-squared differences in the range of 4 eV 2 to 10 4 eV 2 
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when neutrino mixings are of order 10~ 5 < sin 2 29. This result is extremely 
interesting since it covers the range of neutrino masses that give a significant 
contribution to the cosmological energy density. Here the effects of the small 
background neutrino density on the flavor evolution are studied using Eqs. 

It is reasonable to take the supernova to be spherically symmetric. Then 
the density matrix can be parametrized in terms of only 4 quantities: the 
magnitude of the neutrino momentum, k, the radial distance, r, the angle 
between the neutrino velocity and the radial direction, 0, and t. In addition, 
the time dependence in the electron background is negligible since the bulk 
transport of normal matter is nonrelativistic. Thus it is probably reasonable 
to neglect the time dependence everywhere and look for static solutions; i.e. 
p = p(r, k, 0) and p = p(r, k, </>). 

At first glance, numerical simulation would appear to be a straightforward 
way to study the flavor evolution. However as discussed above there are 3 
parameters that must be simultaneously accounted for. In addition, one is 
interested here in what occurs for very small mixing angles — which requires 
very small grid spacings for these three parameters in order to resolve the 
MSW resonances. Both of these considerations make a numerical analysis 
extremely more challenging than anything previously attempted (see e.g. 
[0]). Here a different approach is adopted. 

I start by assuming that the contribution of the neutrino density to the 
potential is negligible to leading order. Then the leading order evolution 
equations are the usual, linear equations that have been well studied (see 
e.g. 0). I then solve these equations and substitute this leading order 
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solution into the full potentials, Eqs. (§-[)]), to get a potential correct to first 
order in the neutrino density. With this corrected potential, the evolution 
equations are again linear and may be easily solved to give the first order 
effects due to a small neutrino background. 

At leading order, expressions for the explicit density matrix solutions can 
be written down. The antineutrinos have little flavor evolution since their 
energy levels are well separated. Assuming the antineutrinos are produced 
at high densities, the general, traceless solution for the antineutrino density 
matrix in the flavor basis is given by 



po(r,k,<f>) = \{fe- fp] 



cos 29 m — sin 29 n 
- sin 29m — cos 29 r , 



(10) 



Here the / Q 's are just the phase space densities as calculated for neutrinos 
that do not mix. For example, for a blackbody spectrum 

fa{k) = [exp (k/T a ) + 1] (H) 

where T a is the temperature of neutrino species v a . Also, 9 m is the leading 
order neutrino mixing angle in matter 

tan29 m (k } r) = = A Sm 26 (12) 

V ; [—2y/2G F N e (r)k + A cos 29] 

and for antineutrinos 9 m is given by the same expression with N e — > — N e . 
There are no explicit phases in Eq. ([HID because they average out in the 
summation of Eq. (^) due to the variation in initial production positions. 
However neutrinos go through a resonance and so have much larger, and 
more complicated, flavor evolution. In the flavor basis the diagonal element 
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of the traceless neutrino density matrix is 

p (r,M) ee = \[fe~ U}{-cos26 m {l-2P c } (13) 



+ 2^P C (1 - P c ) sin 26 m cos(a + (3)} 

and the off diagonal element is 

Po (r,M) eM = ^[/ e -^]{sin20 m [l-2P c ] (14) 
+2y/P c (l- P c )[cos 26 m cos(a + (3) + i sin(a + /?)] } 

For sin 2 26 « 1, the asymptotic^ crossing probability is well described by 
the Landau-Stueckelberg-Zener |13| crossing probability 



P c {k, (b) ~ exp 

where 



71 

" 7 2 



(15) 



7 . (16) 

' 2Pcos2^|ciiV e /Ar e rf s | v ; 

is the adiabaticity parameter and s = k ■ x is the distance along the neutrino 
path. Also a = (s — s )A/2k is the phase generated by vacuum oscillations 
as measured from the resonance point, and (3 is a 'topological' phase acquired 



by going through the resonance [14 



The phases in the leading order neutrino density matrix make it difficult 
to calculate the first order corrections. However the phases vanish and eval- 
uation becomes simple in the adiabatic, P c = 0, and nonadiabatic, P c = 1, 
limits. 



lr The LSZ expression must be amputated above the resonance || but it can reasonably 
be used below the resonance. 
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2.1 The Nonadiabatic limit. 

When the neutrino vacuum mixing is extremely small the flavor evolution 
is extremely small. Then the leading order density matrices are just diago- 
nal in the flavor basis, as is the v — e potential. The effect of the neutrino 
background can then always be incorporated by defining an effective elec- 
tron density. Performing the integration over background neutrinos gives 
the effective neutrino density Nna- Assuming a common neutrinosphere and 
ignoring limb darkening 

N NA (r,<j>) = F(<P,r)[N Ue -N^-N De +N Pt ] (17) 

Here N Ua is the total number density of u a at the neutrinosphere, e.g. for a 
blackbody N Va = 0.0913T^. The geometric factor 

F(4>,r) = -[2-(l + x)cos0][l-x] (18) 

describes the change in the effective neutrino flux outside of the neutri- 



nosphere. Here x = y 1 — (R/r) 2 , and R is the radius of the neutrinosphere 
15] for details on similar geometric factors) .0 



sec 



In the nonadiabatic limit, P{y e —* k, 0) ~ 7- for sin 2 29 << 7 << 1 
and a small neutrino background is fully accounted for by using this expres- 



sion with the replacement N e (r) — ► N e (r) + A^7v^(r, <p) in 7, Eq. (jig) . In 
supernovae the Nna has the same sign as the electron density N e and is at 
most 20%-30% of its size ||. The rate of change of N^a only slightly exceeds 
that of iV e at distances far from the core. Thus the effects of the neutrino- 
antineutrino background on neutrino flavor evolution near the nonadiabatic 
2 When F is averaged over 0, and taking the limit of r >> R, then <F>«|[fl/r] . 
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limit are generally relatively minor. This is apparently the general situation 
the authors of Ref. || envisioned. 



2.2 The Adiabatic limit. 

However, for almost all of the parameter region probed by supernovae nu- 
cleosynthesis, the flavor evolution is extremely adiabatic at leading order in 
the small neutrino background. Substituting the leading order density ma- 
trices with P c = into the potentials and performing the integration over 
background momentum, the diagonal elements of the potential in the flavor 
basis are 

[V^ - K e ] (r, k, 0; po, Po) = ^ cos 26 - V2G F [N e (r) + N Ad (r, 0)] (19) 

where, to leading order in sin 28, 

N Ad {r, 0) ~ F{<t>, r) [-N Ue {q < q r ) + N Ue {q > q r ) + N V/l (q < q r ) 

-N v ^q>q r )+N^-N Pe ] (20) 

Here N Ua (q < q r ) (N Ua (q > q T )) is the number density of v a in equilibrium 
at the neutrinosphere that have momenta less than (greater than) the res- 
onance momentum at radius r, q r = A cos2^/[2v^2G j pA^ e (r)]. The 'reason' 
for the inequalities in NAd is straightforward, the flavor content of the back- 
ground neutrinos depends on whether or not they have gone through their 
resonances. Similarly, the off-diagonal element in the flavor basis is 

K M (r,M;Po,Po) =isin2# [— - V2G F N Ao (r, 0)] (21) 



13 



where, to leading order in sin 26 

iV Ao (r,0) ~ F((f>,r 

-N Ve (q < q r ) + N v (q < q r ) + N p Jq < q r ) - N &e (q < q 



21n|2cot2% r (^| gr -^U (22) 
dq dq 



Here dN a /dq\ qr = qlf a {q r ) / (2vr 2 ) is the derivative of the number density of v a 
in equilibrium at the neutrinosphere, evaluated at the resonance momentum. 
The peak in sin 29 m at the resonance momentum gives the derivative terms, 
while sin 29 m yields the other terms. 

In the adiabatic limit, the effects of the neutrino background on the flavor 
diagonal elements can be incorporated by using in 7 a new, effective electron 
density N' e (r,4>) = N e (r) + NAd(r,4>), analogous to the nonadiabatic case. 
For typical supernovae conditions, NAd is about 2.5 times larger than the 
nonadiabatic neutrino density, Nna, and has the opposite sign from the 
electron density. Thus the effective electron density will be considerably less 
than the physical electron density — especially for large <fr or late times or large 
r. Similarly, the effects of the neutrino background on the flavor off-diagonal 
elements can be incorporated by using in 7 a new, effective vacuum mixing. 
For sin 2 29 « 1 the resonance is narrow so Eq. fl2~T| ) can be evaluated at the 
resonance position, r , to yield 

sin2^(0)^sin2^[l-(^^)] (23) 

The dominant terms in N^o are those with the large logarithmic factor. 
Evaluating at the average muon-neutrino energy, and using the neutrino 
densities given in Ref. § yields that (N Ao /N' e )\ ro > 1 for A < 700 eV 2 , with 
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a maximum of 8 at A ~ 10 eV 2 . Thus the neutrino background is no longer 
'small', for A < 700 eV 2 . 

Of the two limits studied here, the adiabatic limit is the most relevant for 
estimating the effects of a small neutrino background on supernovae nucle- 
osynthesis. This is because when the neutrino background is neglected the 
flavor evolution is extremely adiabatic for almost all of the relevant param- 
eter space, and also because neutrino masses are only probed if significant 
amounts of flavor conversion occur. As A approaches 700 eV 2 from above, 
the 'small' neutrino background increases P c , driving the system away from 
the adiabatic limit. This does not imply that the nonadiabatic limit is then 
reasonable but instead that the oscillatory terms in Eq. fll5l) become impor- 
tant. This is because the MSW effect dramatically enhances the amplitude of 
the oscillatory terms over that of the other flavor off-diagonal terms so that 
only a very small amount of coherence is required for the oscillatory terms 
to dominate. Then the neutrino and antineutrino flavors evolve nonlinearly 
and considerable further study is required in order to reliably describe them. 
For A < 700 eV 2 the connection between r-process nucleosynthesis and cos- 
mologically relevant neutrino masses is enervated. 
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